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Thermal fluctuations strongly modify the large length-scale elastic behavior of crosslinked mem- 
branes, giving rise to scale-dependent elastic moduli. While thermal effects in flat membranes are 
well understood, many natural and artificial microstructures are modeled as thin elastic shells. Shells 
are distinguished from flat membranes by their nonzero curvature, which provides a size-dependent 
coupling between the in-plane stretching modes and the out-of-plane undulations. In addition, a 
shell can support a pressure difference between its interior and exterior. Little is known about the 
effect of thermal fluctuations on the elastic properties of shells. Here, we study the statistical me- 
chanics of shape fluctuations in a pressurized spherical shell using perturbation theory and Monte 
Carlo computer simulations, explicitly including the effects of curvature and an inward pressure. 
We predict novel properties of fluctuating thin shells under point indentations and pressure-induced 
deformations. The contribution due to thermal fluctuations increases with increasing ratio of shell 
radius to thickness, and dominates the response when the product of this ratio and the thermal 
energy becomes large compared to the bending rigidity of the shell. Thermal effects are enhanced 
when a large uniform inward pressure acts on the shell, and diverge as this pressure approaches 
the classical buckling transition of the shell. Our results are relevant for the elasticity and osmotic 
collapse of microcapsules. 
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The elastic theory of thin plates and shells [T], a sub- 
ject over a century old, has recently found new appli- 
cations in understanding the mechanical properties of a 
wide range of natural and artificial structures at micro- 
scopic length scales. The mechanical properties of viral 
capsids [2H1| ) red blood cells [5] , and hollow polymer and 
polyelectrolyte capsules [61-HO] have been measured and 
interpreted in terms of elastic constants of the materi- 
als making up these thin-walled structures. Theoreti- 
cally, models that quantify the deformation energy of a 
two-dimensional membrane have been used to investigate 
the shapes of viral capsids pTlfT5] and their expected re- 
sponse to point forces and pressures [HHIT], as well as 
shape transitions of pollen grains [18] . 

Like its counterparts in other areas of science, such as 
fluid dynamics and the theory of electrical conduction in 
metals, thin shell theory aims to describe the physics of 
slowly varying disturbances in terms of a few macroscopic 
parameters, such as the shear viscosity of incompress- 
ible fluids and the electrical conductivity of metals. De- 
spite such venerable underpinnings as the Navier-Stokes 
equations and Ohm's law, these hydrodynamic theories 
can break down, sometimes in spectacular ways. For ex- 
ample, it is know from mode coupling theory |19j and 
from renormalization group calculations |20| that ther- 
mal fluctuations cause the shear viscosity of incompress- 
ible fluids to diverge logarithmically with system size in 
a two-dimensional incompressible fluid. In the theory of 
electrical conduction, quenched disorder due to impuri- 
ties coupled with interactions between electrons lead to 
a dramatic breakdown of Ohm's law in thin films and 
one-dimensional wires at low temperatures, with a con- 
ductance that depends on the sample dimensions ^TP\. 




FIG. 1. Simulated thermally fluctuating shells, (a) Tri- 
angulated shell with 5530 points separated by average nearest- 
neighbor distance ro with Young's modulus Y — 577e/ro and 
bending rigidity k — 50e at temperature IcbT — 20e, where 
e is the energy scale of the Lennard- Jones potential used to 
generate the disordered mesh, (b) Same as in (a) with ex- 
ternal pressure p = 0.5pc, where Pc is the classical buckling 
pressure. The thermally excited shell has already buckled un- 
der pressure to a shape with a much smaller enclosed volume 
than in (a). 



Even more dramatic breakdowns of linear response 
theory can arise in thin plates and shells. Unlike the 
macroscopic shell structures of interest to civil engi- 
neers, thermal fluctuations can strongly influence struc- 
tures with size of order microns, since the elastic de- 
formation energies of extremely thin membranes (with 
nanoscale thicknesses) can be of the order of the thermal 
energy ksT (where fee is the Boltzmann constant and T 
the temperature) for typical deformations. The statis- 
tical mechanics of fiat solid plates and membranes {i.e. 
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membranes with no curvature in the unstrained state) 
has been studied previously (see [221 El] and references 
therein). Thermal fluctuations lead to scale-dependent 
elastic moduli for flat membranes, causing the in-plane 
elastic moduli to vanish at large length scales while the 
bending rigidity diverges [24l El] . These anomalies arise 
from the the nonlinear couplings between out-of-plane 
deformations (transverse to the plane of the undeformed 
membrane) and the resultant in-plane strains, which are 
second order in the out-of-plane displacements. 

Much less is known about spherical shells subject to 
thermal fluctuations (Fig. la). In fact, the coupling be- 
tween in-plane and out-of-plane modes is significantly dif- 
ferent. Geometry dictates that a closed spherical shell 
cannot be deformed without stretching; as a result, out- 
of-plane deformations provide a first order contribution 
to the in-plane strain tensor [1]. This introduces new 
nonlinear couplings between in-plane and out-of-plane 
deformations, which are forbidden by symmetry in flat 
membranes. We can also consider the buckling of spher- 
ical shells under uniform external pressure, which has no 
simple analogue for plates (Fig. lb). An early exploration 
with computer simulations combined an analysis of the 
elastic energy due to the linear strain contributions of a 
spherical membrane with the nonlinear corrections from 
flat membranes to suggest new scaling behavior for ther- 
mally fluctuating spherical membranes |26j . However, 
an important nonlinear coupling triggered by the curved 
background metric was not considered, nor was the ef- 
fect of an external pressure investigated. Here, we study 
the mechanics of fluctuating spherical shells using per- 
turbation theory and numerical simulations, taking into 
account the nonlinear couplings introduced by curvature 
as well as the effects of a uniform external pressure. 



RESULTS AND DISCUSSION 

Elastic energy of a thin shell 

The elastic energy of a deformed spherical shell of 
radius R is calculated using shallow- shell theory |27j . 
This approach considers a shallow section of the shell, 
small enough so that slopes measured relative to the sec- 
tion base are small. The in-plane displacements of the 
shallow section are parametrized by a two-component 
phonon field Ui(x), i = 1,2; the out-of-plane displace- 
ments are described by a field /(x) in a coordinate system 
x = (2:1,2:2) tangent to the shell at the origin. We focus 
on amorphous shells, with uniform elastic properties, and 
can thus neglect the effect of the 12 inevitable disclina- 
tions associated with crystalline order on the surface of 
a sphere |llj . In the presence of an external pressure p 
acting inward, the elastic energy for small displacements 
in terms of the bending rigidity k and Lame coefficients fi 



and A reads (see Supplementary Information for details) : 
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where the nonlinear strain tensor is 
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Here, d^x = ^Jgdxidx2, where g is the determinant of the 
metric tensor associated with the spherical background 
metric. Within shallow shell theory, 5 ~ 1 (see Supple- 
mentary Information). 

If we represent the normal displacements in the form 
/(x) = /o + /'(x), where /o represents the uniform 
contraction of the sphere in response to the external 
pressure, and /' is the deformation with reference to 
this contracted state so that / d^xf = 0, then the 
energy is quadratic in fields iti, U2 and /o. These 
variables can be eliminated in a functional integral of 
exp(— G'[/', /o, ui, U2]/kBT) by Gaussian integration (see 
Supplementary Information for details). The effective 
free energy Geff which results is the sum of a harmonic 
part Go and an anharmonic part Gi in the remaining 
variable /'(x): 
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where Y = 4/^(/i -I- A)/(2/i + A) is the two-dimensional 



Young modulus and P^^ 



didj/V is the trans- 



verse projection operator. The "mass" term Y{f' /R)"^ 
in the harmonic energy functional reflects the coupling 
between out-of-plane deformation and in-plane stretch- 
ing due to curvature, absent in the harmonic theory 
of flat membranes (plates). The cubic interaction term 
with a coupling constant — y/2i? is also unique to curved 
membranes and is prohibited by symmetry for flat mem- 
branes. These terms are unusual because they have 
system-size-dependent coupling constants. Note that an 
inward pressure [p > 0) acts like a negative i?-dependent 
surface tension in the harmonic term. As required, the 
effective elastic energy of fluctuating flat membranes is 
retrieved for i? -> 00 and p — 0. In the following, we 
exclusively use the fleld /'(x) and thus drop the prime 
without ambiguity. 

When only the harmonic contributions are consid- 
ered, the equipartition result for the thermally generated 
Fourier components fq — J d^x /(x) exp(iq-x) with two- 
dimensional wavevector q are 
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where A is the area of integration in the (2:1,2:2) plane. 
Long-wavelength modes are restricted by the flnite size 
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of the sphere, 



q > 1/R. In contrast to flat 



membranes for which the amphtude of long-wavelength 
[q — > 0) modes diverges as k-BT/{nq'^), the coupling be- 
tween in-plane and out-of-plane deformations of curved 
membranes cuts off fluctuations with wavevectors smaller 
than a characteristic inverse length scale p6] : 



1/4 



7 



1/4 



R 



where we have introduced the dimensionless Fdppl-von 
Kdrmdn number 7 ~ YB? j k We focus here on the 
case 7 1, so ^* <C R- As p approaches pc = 4-\/kF /R^, 
the modes with q = q* become unstable and their am- 
plitude diverges. This corresponds to the well-known 
buckling transition of spherical shells under external pres- 
sure [57]. When p > Pc, the shape of the deformed shell is 
no longer described by small deformations from a sphere, 
and the shallow shell approximation breaks down. 



Anharmonic corrections to elastic moduli 



(See Supplementary Information for details of the calcu- 
lation and the complete dependence on p/pc-) Thus the 
long-wavelength deformations of a thermally fluctuating 
shell are governed by a smaller effective Young's modulus, 
a larger effective bending rigidity, and a nonzero negative 
surface tension even when the external pressure is zero. 
At larger p/pc, however, both the Young's modulus and 
the bending modulus fall compared to their zero tem- 
perature values, and the negative effective surface ten- 
sion determined by p^ gets very large. The complete 
expressions for the effective elastic parameters, includ- 
ing the full p/pc-dependence, show that all corrections 
diverge as p/pc — >■ 1. Furthermore, the effective elas- 
tic constants are not only temperature-dependent, but 
also system size-dependent, since ^/7 oc R. Although the 
corrections are formally small for fceT' <C k, they nev- 
ertheless diverge as i? 00! The thermally generated 
surface tension, strong dependence on external pressure, 
and size dependence of elastic constants are unique to 
spherical membranes, with no analogue in planar mem- 
branes. 



The anharmonic part of the elastic energy, neglected 
in the analysis described above, modifies the fluctuation 
spectrum by coupling Fourier modes at different wavevec- 
tors. Upon rescaling all lengths by £* , it can be shown 
that the size of anharmonic contributions to (|/qP) is set 
by the dimensionless quantities k^T^f^ / k and p/pc- The 
correlation function including the anharmonic terms in 
Eq. [3] is given by the Dyson equation. 



where S(q) is the self-energy, which we evaluate to one- 
loop order using perturbation theory. While (|/qP) can 
be numerically evaluated at any q, an approximate but 
concise description of the fluctuation spectrum is ob- 
tained by expanding the self-energy up to order and 
defining renormalized values Yr, Kr and pn of the Young's 
modulus, bending rigidity and pressure, from the coeffi- 
cients of the expansion: 
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To lowest order in k^T / k and p/pc we obtain the ap- 
proximate expressions (see Supplementary Information 
for details) 
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Simulations of thermally fluctuating shells 

We complement our theoretical calculations with 
Monte Carlo simulations of randomly triangulated spher- 
ical shells with discretized bending and stretching ener- 
gies that translate directly into a macroscopic 2D shear 
modulus Y and a bending ridigity k [28j 2^ ■ (Details are 
provided in Materials and Methods.) Here we study shells 
with 600 < 7 < 35000 and 2 x 10"^ < k^T/n < 0.5. The 
anharmonic effects are negligible at the low end of this 
temperature range. 

The fluctuation spectra of the simulated spherical 
shells are evaluated using an expansion of the radial dis- 
placement field in spherical harmonics [3D]. The radial 
position of a node i at angles (0, 9) can be written as 
rj((/), 6) = Rq + f{<j), 9) with Rq the average radius of the 
fluctuating vesicle. The function /((/>, 9) can be expanded 
in (real) spherical harmonics 
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where Im is the large wavenumber cutoff determined by 
the number of nodes in the lattice {Im + 1)'^ — N j30| . 
The theoretical prediction for the fluctuation spectrum 
including anharmonic effects is {Supplementary Informa- 
tion) 
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Fig. [2] displays our theoretical and simulation results for 
the fluctuation spectrum. At the lowest temperature 
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FIG. 2. Fluctuation spectrum in spherical harmonics. Spherical harmonic amplitude of the shape fluctuations of elastic 
shells plotted against the dimensionless spherical wavenumber I for a shell with R = 40ro,F = 577e/r-o and k — 50e at 
temperatures ksT/ k = 7.4 x 10~* (blue), 0.07 (red) and 0.18 (yellow). The fluctuation amplitudes are scaled by k^T so that 
the spectra at different temperatures would coincide in the harmonic approximation. Each subfigure corresponds to a different 
value of the external pressure: p — (a) and p = 0.2pc (b). The symbols are from Monte Carlo simulations, and the solid lines 



are the theoretical prediction, Eq. 11 using the renormalized elastic constants from perturbation theory (Eqs.lTMl, except for 



the lowest temperature, where the bare elastic constants are used since the anharmonic effects are negligible. 



(corresponding to k-QT ^ j k « 0.1 ^ 1), the spectrum is 
well-described by the bare elastic parameters Y , n and -p. 
At the intermediate temperature (k^T ^ j k k, 10) an- 
harmonic corrections become significant, enhancing the 
fluctuation amplitude for some values of I by about 20%- 
40% compared to the purely harmonic contribution. At 
this temperature, one-loop perturbation theory success- 
fully describes the fluctuation spectrum. However, at the 
highest temperature simulated (k^T^j n, « 24), the an- 
harmonic corrections observed in simulations approach 
50% of the harmonic contribution at zero pressure and 
over 100% for the pressurized shell. With such large cor- 
rections, we expect that higher-order terms in the per- 
turbation expansion contribute significantly to the fluc- 
tuation spectrum and the one-loop result overestimates 
the fluctuation amplitudes. 

Similarly, thermal fluctuations modify the mechanical 
response when a shell is deformed by a deliberate point- 
like indentation. In experiments, such a deformation is 
accomplished using an atomic force microscope [U [S]. 
In our simulations, two harmonic springs are attached 
to the north and south pole of the shell. By changing 
the position of the springs the depth of the indentation 
can be varied (Fig. [3^, inset) . The thermally averaged 
pole-to-pole distance (z) is measured and compared to 
its average value in the absence of a force, (zq). For 
small deformations, the relationship between the force 
applied at each pole and the corresponding change in 
pole-pole distance is spring-like with a spring constant 
k^: {F) = ks{{zo) — {z)). The spring constant is related 
to the amplitude of thermal fluctuations in the normal 
displacement field in the absence of forces by (see Sup- 



plementary Information for the detailed derivation) 
k^T k^T 



2([/(x)]^ 



(12) 



This fluctuation-response relation is used to measure the 
temperature dependence of kg from simulations on fluc- 
tuating shells with no indenters. At finite temperature, 
anharmonic effects computed above make this spring con- 
stant both size- and temperature-dependent: 



k T 

1 -0.0069 — /7 
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Fig. [3^ shows the force-compression relation for a 
shell with R — 20ro and dimensionless temperatures 
kBT^/n = 1.36 X 10"'' and k^iT^/K = 34. The lin- 
ear response near the origin (Fig. |3p) is very well de- 
scribed by fcg measured indirectly from the fluctuations 
in zq at each temperature, Eq. |12| The thermal fluctua- 
tions lead to an appreciable 20% reduction of the spring 
constant for this case. Measuring spring constants over 
a range of temperatures (Fig. [s]:) confirms that the shell 
response softens as the temperature is increased, in agree- 
ment with the perturbation theory prediction. We note, 
however, a small but systematic shift due to the finite 
mesh size of the shells, an approximately 5% effect for 
the largest systems simulated here. At the higher tem- 
peratures (k-QT^jn > 20), the measured spring con- 
stants deviate from the perturbation theory prediction, 
once again we believe due to the effect of higher-order 
terms. 

We also simulate the buckling of thermally excited 
shells under external pressure. When the external pres- 
sure increases beyond a certain value (which we identify 
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FIG. 3. Temperature dependence of response to point forces, (a) Force-compression curves for simulations of indented 
shells (symbols) with R = 20ro, Y = 577e/ro and k = 50e at low [ksT/K = 2 x 10~^) and high [ksT / n = 0.5) temperature. 
The lines show the expected linear response at small deformations with the spring constant fcs measured independently from 
fluctuations in zo (fcs — 29.15e/ro for ksT/K = 2 x 10~^, fcs = 23.63e/r-Q for k^T/n — 0.5). For indentation depths larger 
than 1 — {z}/{zo} ~ 0.05, the regions around the poles become inverted and the response becomes nonlinear. Inset: schematic 
showing the definition of zq (the pole-to-pole distance in the absence of indentations) and z (pole-to-pole distance following an 
indentation imposed by harmonic springs whose free ends are brought close together) for a snapshot of the fluctuating shell, (b) 
Blow-up of the boxed region near the origin in (a), highlighting the linear response regime, (c) Spring constants extracted from 
fluctuations for shells with three different radii as a function of temperature, rescaled by the classical result for linear response 
of thin shells at zero temperature. The dashed line shows the perturbation theory prediction, Eq. |13[ The low-temperature 
spring constant deviates from the classical result due to a finite mesh size effect which falls with increasing R (increasing mesh 
size). 



as the renormalized buckling pressure) , the sheU coUapses 
from a primarily spherical shape (Fig. la) to a shape with 
one or more large volume-reducing inversions (Fig. lb). 
For zero temperature shells, this buckling is associated 
with the appearance of an unstable deformation mode in 
the fluctuation spectrum. At finite temperature, the ap- 
pearance of a mode with energy of order k^T is sufficient 
to drive buckling. Anharmonic contributions, strongly 
enhanced by an external pressure, also reduce the effec- 
tive energy associated with modes in the vicinity of q* 
primarily due to the enhanced negative effective surface 
tension Pu.R/'2. (see Eq. [s]). As a result, unstable modes 
arise at lower pressures and we expect thermally fluc- 
tuating shells to collapse at pressures below the classi- 
cal buckling pressure Pc- This is confirmed by simula- 
tions of pressurized shells (Fig. |4]). When anharmonic 
contributions are negligible [k-oT^/K <C 1), the buck- 
ling pressure observed in simulations is only ~ 80% of 
the theoretical value because the buckling transition is 



highly sensitive to the disorder introduced by the ran- 
dom mesh. Relative to this low temperature value, the 
buckling pressure is reduced significantly when k-^T^ / n 
becomes large. 



Conclusion and outlook 

In summary, we have demonstrated that thermal cor- 
rections to the elastic response become significant when 
kBTyf^/K ^ 1 and that first-order corrections in k^T/n 
already become inaccurate when k^T^ / k > 20. Hu- 
man red blood cell (RBC) membranes are known exam- 
ples of curved solid structures that are soft enough to 
exhibit thermal fluctuations. Typical measured values 
of the shear and bulk moduli of RBC membranes cor- 
respond to F « 25 ^N/m [3131], while reported values 
of the bending rigidity k, vary widely from 6 k^T to 40 
k^T [SllSlj. Using an effective radius of curvature i? w 7 
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FIG. 4. Temperature dependence of the buckling pres- 
sure. Buckling pressure for simulated shells at various radii 
and temperatures, normalized by the classical {i.e. zero tem- 
perature) critical buckling pressure Pc for perfectly uniform, 
zero temperature shells with the same parameters. For all 
shells, Ytq/k = 11.54. In separate sets of symbols, we either 
vary the shell radius over the range 7.5 < R/ro < 55 while 
keeping the temperature constant {ksT = 2 x 10~^k, blue 
circles; k^T = 0.4k, yellow squares) or vary the temperature 
over the range 2 x 10~* < ksT/K < 0.4 while keeping the 
radius constant at _R = 20ro (red triangles). The parame- 
ter ksTy/^/K sets the strength of anharmonic corrections for 
thermally fluctuating shells. The inset shows the l/R^ de- 
pendence of the buckling pressure as the radius is varied, for 
shells at low and high temperature. 



fim [5] gives k-gT^/K in the range 2-35. Thus, RBCs 
could be good candidates to observe our predicted ther- 
mal effects, provided their bending rigidity is in the lower 
range of the reported values. 

For continuum shells fabricated from an elastic mate- 
rial with a 3D Young's modulus E, thickness h and typ- 
ical Poisson ratio w 0.3, kBT^/n w 100RkBT/{Eh'^). 
Hence very thin shells with a sufficiently high radius- 
to-thickness ratio [R/h) must display significant ther- 
mal effects. Polyelectrolyte [SI and protein-based [55] 
shells with R/h k, 10"^ have been fabricated, but typi- 
cal solid shells have a bending rigidity k, several orders 
of magnitude higher than k^T unless /i ^ 5 nm. Micro- 
capsules of 6 nm thickness fabricated from reconstituted 
spider silk with i? sa 30 /im and E k, \ GPa have 
k-^T^^ / K « 3, and could exhibit measurable anharmonic 
effects. 

Thermal effects are particularly pronounced under fi- 
nite external pressure — an indentation experiment car- 
ried out at p = Pc/2 on the aforementioned spider silk 
capsules would show corrections of 10% from the classical 
zero-temperature theory. For similar capsules with half 



the thickness, perturbative corrections at p = Pc/2 are 
larger than 100%, refiecting a drastic breakdown of shell 
theory because of thermal fiuctuations. The breakdown 
of classical shell theory explored here points to the need 
for a renormalization analysis, similar to that carried out 
already for flat plates |34j . 

MATERIALS AND METHODS 

Monte Carlo Simulations of randomly triangulated 
shells 

A random triangulation of radius i?o is constructed by 
distributing N nodes on the surface of a sphere with the 
required radius. The first two of these nodes are fixed at 
the north and the south pole of the sphere whereas the 
positions of the remaining N ~ 2 nodes are randomized 
and equilibrated in a Monte Carlo simulation. During 
this equilibration process the nodes interact via a steeply 
repulsive potential (the repulsive part of a Lennard Jones 
potential). After equilibration, when the energy has 
reached a constant value on average, the simulation is 
stopped and the final configuration is 'frozen'. The neigh- 
bours of all nodes are determined using a Delaunay tri- 
angulation [55] . The spherical configurations as well as 
the connection lists are used in further simulations. 

In subsequent simulations nearest neighbours are per- 
manently linked by a harmonic potential giving rise to a 
total stretching energy |36j, 

^^. = ^E(I''^^--'^°I')' (A.14) 

where the sum runs over all pairs of nearest neighbours, 
rij is the distance between two neighbours and r^j the 
equilibrium length of a spring. The equilibrium length 
rfj is determined at the start of the simulation, when 
the shell is still perfectly spherical and thus the stretch- 
ing energy vanishes for the spherical shape. The spring 
constant k is related to the two-dimensional Lame coef- 
ficients \ = n = a/3/c/4 and the two-dimensional Young 
modulus Y = 2fc/\/3 [55]. 

The mean curvature (more precisely, twice the mean 
curvature) at node i is discretized using [501 [13 HI] 

= -n,.E^(r,-r,) (A.15) 

where is the surface (unit) normal at node i (the 
average normal of the faces surrounding node i), di — 
'^ij^ij is the area of the dual cell of node i, Oij = 
lij[cot9i + cot02]/2 is the length of a bond in the dual 
lattice and lij — jr; — r^l is the distance between the 
nodes i and j. The total curvature energy is, 

Sb- ^E'T,(if,-i/o)' (A.16) 
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with K the bending rigidity and Hq the spontaneous cur- 
vature at node i. In all simulations Hq = 2/Ro (since Hi 
is twice the mean curvature) . In the cases of elastic shells 
under pressure a term PV is added to the Hamiltonian 
where P is the external pressure and V the volume of the 
shell. 

Similar elastic networks with stretching and bending 
potentials have been studied in relation to the stability of 
membranes, icosahedral and spherical shells that contain 
defects [HI [HI EH IMl IMl SO] or defect scars [38l M- 
H5] as well as for the deformation of icosahedral viruses 
[TIHTO] and the crumpling of elastic sheets [55]. 

Simulations are performed for shells of 5530 {Rq — 
20 ro), 22117 (i?o = 40 ro) and 41816 {R = 55 tq) 
nodes. The Hookean spring constant and the bend- 
ing rigidity are taken such that the shells have Foppl- 
von-Karman numbers in the range 650 < 7 < 35000 
and that the dimensionless temperature is in the range 
2 X 10"'^ < ksT/n < 0.5. Monte Carlo production runs 
consist typically of 1.25 x 10^ Monte Carlo steps where in 
a single Monte Carlo step an attempt is made to update 
the positions of all nodes once on average. Configurations 



were stored for analysis typically every Ng-, 



2000 



Monte Carlo steps. For the largest system (41816 nodes), 
such a run took about 700 days of net CPU time spread 
over several simultaneous runs in a Linux cluster of Intel 
XEON X5355 CPUs. For the smaller shells, the compu- 
tational time scaled down roughly linearly with system 
size. 



The fluctuation spectrum from computer simulations 



shell. This leads to an additional term in the Hamiltonian 
Vs = h {zf - z^f /2 + h {zf - z^f /2 where h = t^/rl 
is the spring constant of the indenter. Here, one end of 
the springs, at positions z^ and z^, is attached to the 
vertices at the north and south pole, respectively. The 
positions of the other end of the springs, at zf and zf , 
are fixed externally and determine the indentation force 
and depth, as indicated in Fig. S4. 

By changing zf and zf , the depth of the indentation 
can be varied. After the springs are fixed a certain dis- 
tance apart, the thermally average pole-to-pole distance 
(z) is measured and compared to its value in the absence 
of a force, (zq). The instantaneous force at the poles is 
calculated from the instantaneous extension of the har- 
monic springs after each iVgamp Monte Carlo steps; ther- 
mal averaging then determines the average corresponding 
to (z). This provides the force-indentation curves in Fig. 
3(a-b). 

It is very difhcult to unambiguously identify the linear 
regime in the force-indentation curves. Extracting the 
effective spring constant of shell deformation kg from a 
linear fit in the small indentation region is subject to inac- 
curacies and sensitivity to the number of points included 
in fitting. Instead, we extract the spring constants of 
thermally fluctuating shells by using a relation between 
ks and the fluctuations in zg (see Supplementary Infor- 
mation for derivation): 



(A.17) 



For a particular configuration of a simulated shell, the 
coefficients Aim of the expansion of the radial displace- 
ments in spherical harmonics (Eq. 10 1 are determined 
by a least squares fit of the node positions to a finite 
number Im of (real) spherical harmonics. In practice we 
have used Im — 26 as the upper wavenumber cutoff for 
all simulations. At each temperature and pressure, this 
procedure is repeated for about 10000 independent con- 
figurations and the results averaged to obtain the curves 
presented in Fig. [2] 



Simulations of shells indented by point-like forces 

To perform indentation simulations, two harmonic 
springs are attached to the north and south pole of the 



This procedure was used to measure the temperature- 
dependent spring constants in Fig. 3c. 
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FIELDS AND STRAINS IN SHALLOW SHELL THEORY 



Z 

A 




FIG. SI. The coordinate system in shallow shell theory. A section of the undeformed sphere is shown with the (x\,x-2) plane 
tangential to it at the origin. The red arrows show the directions into which displacements Mi(a::i, 3:2), U2(x\,x-2) and f{x-i_,X2) 
are decomposed at a particular point in the coordinate plane. 



We describe the deformations of the sphere using shallow shell theory which we summarize here. We follow the 
presentation by Koiter and van der Heijden [T]. A shallow section of the sphere is isolated and Cartesian coordinates 
{xi,X2) are set up to define a plane that just touches the undeformed sphere at the origin and lies tangent to it; 
the z axis is thus normal to the sphere at the origin (Fig. SI). We use the Monge representation to parametrize the 
undeformed shell by its height z = Z{xi,X2) above the plane, where Z{xi,X2) is the undeformed state corresponding 
to a sphere of radius R with its center located on the z-axis above the {xi, X2) plane; 

^(-l--2)=i?(^l-yi-§-§) (SI) 

The assumption in shallow shell theory is that the section of the shell under consideration is small enough that slopes 
diZ ^ xi/R and 822 ^ X2/R measured relative to the {xi,X2) plane are small. (Partial derivatives are denoted by 
d/dxi = di.) Then the undeformed state is approximately parabolic in xi and X2, 

Z{x,,X2)^^^^. (S2) 

Deformations from this initial state are quantified via a local normal displacement f{xi,X2) perpendicular to the 
undeformed surface and tangential displacements ui{xi,X2) and U2{xi,X2) within the shell along the projections of 
the xi and X2 axes on the sphere respectively. In terms of these fields, a point {xi,X2, Z{xi,X2)) in the undeformed 
state moves to (a^i + ui — fdiZ, X2 + U2 — fd2Z, Z + f ) to lowest order in the slopes diZ — Xi/R. The strain tensor 
is defined by the relation between the length ds' of a line element in the deformed state and the corresponding line 
element length ds in the undeformed state [5]: 

{ds')"^ = ds^ + 2uijdxidxj . (S3) 

With this definition and neglecting terms of order (diZ)^ and their derivatives, we find the nonlinear strain tensor 
used in the main text, 

1 / 
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The stretching energy is then given by [2] 

Gs^lJdS [2^4 + AuL] , (S5) 

where /i and A are the Lame coefficients and dS is an area element. 

We also include a bending energy of the Helfrich form [5] that penalizes changes in local curvature: 

Gb^^l dSiH-Hof, (S6) 

where k is the bending rigidity, H the mean curvature and Hq the spontaneous mean curvature (which we take to be 
equal everywhere to the curvature 2/R of the undeformed shell). For a shallow section of the shell, the local curvature 
can be written in terms of the height field Z{xi,X2) + f(xi, X2) as 

H = W^{Z + f) = ^+V'f, (S7) 

where = dn + 822 is the Laplacian in the tangential coordinate system. Finally the energy due to an external 
pressure p equals the work done. 



W = -pJ dSf. (S8) 

The area element is dS = dxidx2/ \/^ — {xf + x^/i?^ ~ dxidx2 when terms of order [xi/R)"^ and above are 
neglected. Summing the stretching, bending and pressure energies leads to the elastic energy expression G = Gs + 
Gh + W, Eq. (1) in the main text. 

Since we are restricted to a shallow section of the shell, the theory is strictly applicable only to deformations whose 
length scale is small compared to the radius R. The typical length scale £ of deformations can be obtained by balancing 
the bending and stretching energies Gb and Gg discussed above. Upon noting that the stretching free energy density 
in a region of size £ is Gs ^ Y{f / R)'^ , where F is a typical elastic constant, and Gb ~ k/^/^^: we recover the Foppl-von 
Karman length scale introduced in the main text, 

^ (S9) 



where the Foppl-von Karman number is 7 = YR^ / k. More sophisticated calculations (sketched below) show that the 
relevant elastic constant is the 2D Young's modulus, Y = 4:fi{iJ. + A)/(2/i + A). 

For a shell made up of an elastic material of thickness h, taking Y and k from the 3D Young's modulus of an 
isotropic solid within thin shell theory provides the estimate 7 « 10{R/h)'^ [51. For shallow shell theory to be valid, 
we need £* <^ R. Hence shallow shell theory is valid when 7^1 i.e. R ^ h, which is precisely the limit of large, 
thin curved shells which are most susceptible to thermal fluctuations. This agreement between shallow shell theory 
and more general shell theories that are applicable over entire spherical shells has been discussed by Koiter [4 in the 
context of the response of a shell to a point force at its poles. Shallow shell theory was also used to study the stability 
of pressurized shells by Hutchinson [5]. In both cases, shallow shell theory was shown to be valid for thin shells such 
that h/R ^ 1. Since thermal fluctuations are only relevant for shells with radii several orders of magnitude larger 
than their thickness, shallow shell theory is an excellent starting point for the extremely thin shells of interest to us 
here. 



ELIMINATION OF IN-PLANE PHONON MODES AND UNIFORM SPHERICAL CONTRACTION BY 

GAUSSIAN INTEGRATION 

A spherical shell under the action of a uniform external pressure that is lower than the critical buckling threshold 
responds by contracting uniformly by an amount /q. The out-of-plane deformation field can then be written as a sum 
of its uniform and non-uniform parts, 

/(x) = fo + /'(x) = fo + Yl /qe"^''■^ (SIO) 

where /'(x) represents the contribution to the field from its q 7^ Fourier components. (In this section, for ease 
of presentation we use the normalization /q = ^ J d^x f{x)e^'^'^, where A is the area of integration in the (xi,X2) 
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plane. The inverse transform is then /(x) = X]q /q^"^'' With this decomposition, / d?xf'{x) = and thus only 
fo contributes to the pressure work W. On the other hand, only /' contributes to the nonlinear part of the strain 
tensor. Hence the elastic energy G — Gb + Gs + W defined above is harmonic in the in-plane phonon fields ui(x) and 
M2(x) as well as the uniform contraction /q. To analyze the effects of anharmonicity, it is useful to eliminate these 
fields and define an effective free energy [6], 



Gctf[f] = -kBT\nl^J'Du{xi,X2) J dfoe~ 



•G[/',/o,«l,«2]/fcBT] 



(Sll) 



To carry out the functional integrals in Eq. (Sll) for a fixed out-of-plane displacement field /'(x), the strain tensor 
Uij must also be separated into its q = and q ^ components: 



q#0 



i f 

- {q^Uj{q) + qjUi{q)) + A,j{q) - 5,j-^ 



-?q-x 



(S12) 



where 



(S13) 



The uniform part of the strain tensor has the following components: 



7/" 
"12 



n2 



^ii(O)- 



/o 



i?' 

^22(0) 
^12(0). 



(S14) 



Here, u^^- are the uniform in-plane strains that are independent of /q. This restriction implies that + — 
because a simultaneous uniform in-plane strain of the same sign in the xi and X2 directions corresponds to a change 
in radius of the sphere and thus cannot be decoupled from /q. Hence in addition to /o and U12, there is only one more 
independent degree of freedom, Am° = ufi — U22, that determines the uniform contribution to the strain tensor. 

Finally we perform the functional integration in Eq. (Sll I over the phonon fields Ui as well as the three independent 
contributions to the uniform part of the strain tensor — /o, U12 ^^'^ Au'^. The resulting effective free energy is, upon 
suppressing an additive constant. 



Goff 



d^x 



2 i-/\2 



- (^P^d f'd f- — 
2 I 2^'^ R 



[An(0) 



^22(0)] 



(S15) 



where = Sij — didj/V"^ is the transverse projection operator. Note that as a result of the integration the Lame 
coefficients /i and A enter only through the 2D Young's modulus Y — 4:fi{ijL + A)/(2/i -I- A). Finally, substituting 



All(0)+A22(0) 



1 

2A 



d'x [idjr+{d2f'r 



1 

2A 



d^x\Vf 



l\2 



(S16) 



in Eq. ( S15 1 gives the effective free energy used in the analysis, Eq. (3) in the main text. In the following, we drop the 



prime on the out-of-plane displacement field since /q has now been eliminated. When only the harmonic contributions 
are considered, the equipartition result for the thermally generated Fourier components /q = J d^a; /(x) exp(zq • x) 
with two-dimensional wavevector q are 



(/q/q')o 



2 9 + fl2 



(S17) 



where A is the area of integration in the {xi,X2) plane. This harmonic spectrum [Eq. (4) in the main text] takes on 
corrections due to the anharmonic terms that are calculated in the next section. 
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'^1 \, ./ '^3 



1^2 ^4 



^[P^j(kl + k2)fcHfc2,][P;™(k3 + U)k3lk4m]{27rfd^]^i + ka + kg + k4) 



k. 



1 



k2 



2R 



'3 — [f^](ki + k2)kuk2jK2TrfS\ki + ka + kg) 



FIG. S2. The bare propagator for /(q) and the vertices arising from the non-quadratic terms in GcH- The slashes on specific 
legs denote spatial derivatives. Pij(q) = Sij — qiqj/q^ is the transverse projection operator in momentum space. Note an 
unusual feature of this graphical perturbation theory: the system size, i.e. the sphere radius R, enters explicitly both in the 
propagator and as a coupling constant in the third order interaction vertex. 






FIG. S3. One-loop corrections to the two-point height-height correlation function in momentum space. Note that in calculating 
the self-energy, the external propagators are not included; i.e. they are "amputated". The contribution in a resembles that for 
membranes with a flat ground state, except for the ii-dependent pressure and mass terms in the propagator. The nonlinear 
corrections in b, however, arise from a cubic coupling constant proportional to 1/R and are unique to the spherical geometry. 

ONE-LOOP CONTRIBUTIONS TO THE SELF-ENERGY 

Here we describe the self-energy used to calculate the leading anharmonic corrections to the fluctuation spectrum 
in the main text. The Feynman rules obtained from the effective free energy Gcff [/] are summarized in Fig. S2. 
Henceforth, Fourier components are deflned as in the main text: /q = / d^x /(x) exp(iq • x) with two-dimensional 
wavevector q. The inverse Fourier transformation of the out-of-plane deformation field is 



A ^ 

where A is the area of integration in the {xi,X2) plane and the sum is over all allowed Fourier modes. The one- 
loop contribution to the self-energy E(q) due to the anharmonic three-point vertex (cubic term in the energy) and 
the four-point vertex (quartic term) are summarized in Fig. S3. Fig. S3a is also present in the calculation for flat 
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membranes [5], and provides a contribution 



(2^)2 «|q + k|4-£f|q + k|2 + ^ 

to the self-energy. Fig. S3b involves two-vertex terms arising from the cubic coupling unique to shells with curvature 
(note that, despite "amputation" of the propagator legs, the diagrams are distinct because the slashes decide the 
momentum terms which survive various index contractions in addition to determining the momentum of the transverse 
projection operator introduced at each vertex). The net contribution to the self-energy from the four diagrams in 
Fig. S3b is 



y2 



^|q + k|4-2^|q- 



k|2 + y\ (^.k^ - vRq^ 



Y 



(S20) 



-[P^^{ci)hk,Y + [P;^(k)g,g,f + [P]^^{\^)q,q.j\[P^,^{\^ + Ci)qiq„;\ + 2[P'^^{\^)q,q,][P^,^{ci)kik„ 



While the inverse of the harmonic correlation function, Eq. (S17), only contains terms of order q^, q^, the 



one-loop corrections to the spectrum [Eqs. (S19 S20)] generate terms with these powers of q as well as terms of order 
q^ and above in the full inverse fluctuation spectrum. If we keep only terms of order and below in the calculation 
of the one-loop fluctuation spectrum, we can provide an approximate description of the low-q behaviour of the shell 
in terms of effective elastic constants: 



AfcBT(|/q^0|'> 



PkR 2 

-^<1 



i?2 



(S21) 



where Y^, Kr and are the effective Young's modulus, bending rigidity and dimensionless pressure respectively. At 
long length scales, probes of the elastic properties of thermally fluctating shells would provide information of these 
effective elastic constants rather tha n the "ba re" constants Y, k and p that describe the zero-temperature shell. Upon 
expanding the integrands in Eqs. (S19 S20l to 0{q'^) the momentum integrals can be carried out analytically to 
obtain: 



3 fceT 



1287r K (1-772)3/2 



^7/-\/l — -|- TT — COS ^ T]j 



(S22) 



1 ksT ^ 



(-1699 + 3758r;2 - 21047?^ 



307207r K (1-772)7/2 

+ 15(61 - 288r/2 + 41677^ - 19277^) (tt - cos"^ 77) 



(S23) 



Vn. = V + 



1 knT 



15367r K (1 _ 7y2) 



2^5/2 



v/1-772 (64 - 67772) + 3 (2I77 - 2277^) (tt - cos^^ : 



(S24) 



where we have defined a dimensionless pressure 77 = p/pc and p^ = 4%/ kY /R^ is the classical buckling pressure of 
the shell. We see explicitly that the quantities diverge in the limit 77 — > 1. To lowest order in the external pressure, 
we have 



Y^^Y 



Pr~P 



1 - 



3 keT 



256 K 
247r K 



^ 1 + 



4 p 



and 



61 fceT 

1 + 1 

4096 K ^ ' ^ 



np, 

637r p 
128pc 

1568 p 
9157r pc 



(S25) 
(S26) 

(S27) 
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These are the approximate renormalized elastic quantities tabulated in Eqs. (7-9) in the main text. 

In evaluating the above expressions, the momentum integrals in Eqs. (S19 S20l must strictly speaking be carried 
out over the phase space of all allowed Fourier modes /(k) of the system, which go from some low-fc cutoff k^i^ ^ ^/R 
to a high-fc cutoff set by the microscopic lattice constant. However, since all integrals converge in the ultraviolet limit 
fc — > oo, the upper limit of the fc-integrals can be extended to oo. The integrals are well-behaved at low momenta 
due to the mass term ^ Y/R^ in the propagator. Hence we carry out the momentum integrals over the entire 
two-dimensional plane of k. The excess contribution to the self energy by including spurious Fourier modes with 
< A: < 1/i?, i.e. for wavevectors less than the natural infrared cutoff k^in ^ 1/i?, gives rise to an error of roughly 
1/^/7 which is negligible for extremely thin shells. This correction is of similar magnitude to the errors introduced 
by using shallow shell theory (which is inaccurate for the longest-wavelength modes with wavevector k ^ ^/ R) which 
are also negligible in the thin-shell limit. 



CALCULATION OF FLUCTUATION SPECTRUM WITH SPHERICAL HARMONICS 



While the perturbation theory calculations were carried out using a basis of Fourier modes in a shallow section 
of the shell to decompose the radial displacement field, the fluctuation spectrum is most efficiently measured in 
simulations using a spherical harmonics expansion. To compare the simulation results to the expected corrections 
from perturbation theory, we use the description of the shell in terms of the effective elastic constants Yr, Kb. and p^, 
Eqs. (7-9) in the main text. 

Consider a spherical shell of radius R with bending rigidity k and Lame coefficients A and fi, experiencing a tangential 
displacement field u = (ux-, Uy) and a radial displacement field /. Like any smooth vector field, u can be decomposed 
into an irrotational (curl-free) part and a solenoidal (divergence-free) part: u = V^* + v, where the scalar function 
^ generates the irrotational component and v is the solenoidal component. Upon expanding / = ^ AimRYj^ and 



^ = Si m BimR^Yi"^ in terms of spherical harmonics Yf 
order in the fields is given by [7] 



), the elastic energy of the deformation to quadratic 



G = R' 



I, in 



2 R^ 



2K 



Im 



+ -l{l + 1) [{K + (1)1(1 + 1) - 2/i] B'L } + G,oi(v), 



(S28) 



where K = \ + fi is the bulk modulus. The solenoidal component v does not couple to the radial displacement field 
and provides an independent contribution Gsoi which is purely quadratic in the field v. 

To this elastic energy, we also add the surface energy-like contribution Gs = —{pR/2)AA due to the "negative 
surface tension" —pR/2 present in the shell when it is uniformly compressed in response to an external pressure p. 
Here AA is the excess area due to deformations about the average radius. In terms of spherical harmonic coefficients, 
this area change can be written j8] 



AA^R^ A 



i>l,m 



1 



^(/ + 1) 



(S29) 



As we did for the elastic energy in shallow shell theory, we can now integrate out the quadratic fluctuating quantities 
Bijn and the solenoidal fleld v to obtain an effective free energy in terms of the radial displacements alone: 



GcS — 



E 

l>l,m 



n{l + 2f{l-lf 



~pR 



1 



^(/ + 1) 
2 



Aii{(i + \){P + 1-2) 
{2fi + X){P + I) - 2fi 



42 



(S30) 



The fluctuation amplitude is obtained via the equipartition theorem: 

1(1 + 1) 



2\-l 



kBT{\Aim\% 



k{1 + 2)^(1 
k{1 + 2)^{1 



1 



1) 



pR^ 
pR^ 



1 



1 



2 

l{l + l) 



4M(/i + A)(/2 + / + 2)^2 



{2fi + X){P + I) - 2(1 
^ -R^ 



1 



Y 



(S31) 



where Y = ifidj, + A)/(2/i + A) is the 2D Young's modulus introduced earlier. The effect of anharmonic contributions 
to the fluctuation spectrum can now be calculated by using the effective temperature-dependent quantities Kr, Kr 
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and Pr in place of the bare elastic constants in the above expression. However, the last term in Eq. (S31| also 



requires knowledge of the thermal corrections to the Lame coefficient fi which was eliminated in the shallow shell 
calculation when the tangential displacement fields were integrated out. For the discretized stretching energy used in 
the simulations, we have /i = 35^/8. If we assume that this relationship is not significantly changed by the anharmonic 
corrections to one- loop order, the n /Xr « 3Ye/8. Upon substituting this approximation together with the other 
effective elastic parameters in Eq. (S31), we find 



kBT{\Ai„fr^ + 2f{l - If - p^R^ 

which is the same as as Eq. (11) in the main text [9]. 



1 



/(/ + 1) 



3{P + l-2) 



3{P + l)-2 



(S32) 



LINEAR RESPONSE OF THE SHELL TO POINT FORCES 



We calculate the response of the shallow shell to a point force at the origin, corresponding to a force field h{x) 
FS'^{x.). The Fourier decomposition of this force field is 



F, for allq. 



(S33) 



The linear response of the deformation field / to this force is related to its fiuctuation amplitudes in the absence of 
the force, (|/qP)h=oi by the fluctuation-response theorem: 



iU) 



(l/q 



AkBT ''^ 



(l/ql 



)h=0 



AkeT 



F. 



The inward deflection at the origin is then 

(/(x = 0)) = i^(/q) = 



F 



A^kgT 



(S34) 



(S35) 



This can be related to (/^), the mean square fluctuations of the deformation field in real space which is a position- 
independent quantity in the absence of nonuniform external forces: 



([/(x)]'>.=o 



q q' 



(S36) 



From Eqs. (S35) and (S36), we obtain 



(/(x = 0)) = — cf) 
knT 



(S37) 



This equation relates the depth of the indentation due to a force F at the origin to the mean square fluctuations of 

the deformation field / in the absence of such a force. 

When only harmonic contributions are considered, Eq. (S17) gives us the mean square amplitude (|/qP)o = 
AksT / {nq^ — pRq^ /2 + Y / R^) in terms of the elastic constants and external pressure. Upon taking the continuum 
limit of the sum over wavevectors — > ^ / (P'qji^'Kf , we can calculate the fluctuation amplitudes exactly: 



knT 



(27r)2 



nq-^ 



2 9 + 



RkeT 1 + I sin~^ 



(S38) 



where r\ = p/pc = pi?^/(4V kY') is the dimensionless pressure, and ?7 < 1, i.e. we restrict ourselves to pressures below 
the classical buckling pressure. From Eqs. (S37) and (S38l, we get the linear relation between the indentation force 
and the depth of the resulting deformation: 



F = 



r^(/(x = 0)). 



- Sm 77 



(S39) 
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The temperature drops out and we obtain a result valid for T = shells as well. The expression reproduces the 
well-known Reissner solution |10j for the linear response of a spherical shell to a point force when 77 = 0, and also 
reproduces the recent result from Vella et al [TT] for indentations on spherical shells with an internal pressure when 
77 < 1. At finite temperatures, however, anharmonic effects contribute terms of order (ksT)^ and higher to (/^), 
making the response temperature-dependent. 

In the simulations, the shells contract by a small amount due to thermal fluctuations, even in the absence of external 
forces. Thus, indentations are measured relative to the thermally averaged pole-to-pole distance of the shell at finite 
temperature, (zq) < 2R. Equal and opposite inward forces are applied to the north and south poles of the shell to 
maintain a force balance (see details in the Materials and Methods section of the main text) and the resulting average 
pole-to-pole distance, (z), is measure d. T his corresponds to an average indentation depth of {{zq) — {z))/2 at each 
pole, with associated force [from Eq. (S37)] 



keT ((zq) - (z)) _ 
F = = h((zo} - {z}}, 



iP) 2 

i.e. the shell as a whole acts as a spring with spring constant 



2(/2 



At r = 0, we have 



ka — 



R l + fsin"S' 



(S40) 



(S41) 



(S42) 



in particular, fcg — 4%/ kY j R in the absence of external pressure. Anharmonic contributions change the fluctuation 
amplitude (/^) and hence the linear response. To lowest order in temperature, the effects of anharmonic contributions 



can be obtained by using the renormalized elastic constants calculated using perturbation theory [Eqs. (S25)-(S27|] 
in Eq. ( S42 ) and keeping terms to 0(T) . In particular, even if the bare pressure p = 0, the renormalized dimensionless 



pressure is nonzero and affects the spring constant, as do the temperature-dependent effective elastic moduli. The 
result in this case is 



fc,(r>o) 



k T 

1 - 0.0069^^^7 



(S43) 



This is the theoretical prediction quoted as Eq. (13) in the main text. 



MEASURING THE EFFECTIVE SPRING CONSTANT FROM MONTE CARLO SIMULATIONS 



We extract the spring constants of thermally fluctuating shells for Fig. 3c in the main text by using the relation 



between kg and fluctuations in the transverse displacement field / [Eq. (S41)]. It is straightforward to measure the 
average pole-to-pole distance of the fluctuating shell in the absence of external forces, (zq) = (i? — /n ^ /s)i where /n 
and /s are the inward displacements at the north and south poles respectively. Since the displacements at the poles 
are expected to be independent of each other, the mean squared fluctuations in zq are closely related to the mean 
square fluctuations in /: 

{zl)-{z,)^^2{f). (S44) 
The spring constant can thus be measured indirectly from the fluctuations in the pole-to-pole distance using Eq. ( [S4l] ): 



2{P) {zD - {z,Y 



(S45) 
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FIG. S4. Illustration of the indentation simulation of a randomly triangulated shell. Shown here are vertical cuts through a 
shell of radius Ro = 20 ro at dimensionless temperature kBT/n^ = lO"'' (a & b) and ksTjK^ = 15 (c). A shell at low 
temperature (a) is indented by two harmonic springs (dark blue lines) attached to the north («^) and south (z^) pole of the 
shell. Fixing the springs at a separation Zi — zf — zf leads to a pole separation z. The green points indicate the positions zf 
and zf of the fixed end points of the springs, (b) A close-up of the north pole of the shell displayed in a. The configuration 
contains a minority of 5 and 7-fold coordinated vertices in addition to 6-fold coordinated ones, (c) Illustration of a fluctuating 
shell at Ro = 20 ro and ksT/K^ = lO"'^. 
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